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Abstract 

A constructive approach to get the reduced row echelon form of a given 
matrix A is presented. It has been shown that after the kth step of the 
Gauss- Jordan procedure, each entry a k j(i ^ j,j > k) in the new matrix 
A k can always be expressed as a ratio of two determinants whose entries 
are from the original matrix A. The new method also gives a more general 
generalization of Cramer's rule than existing methods. 

Key words: Gauss- Jordan Elimination, Cramer's rule, Determinants 



1. Introduction 

Gauss- Jordan elimination is a variation of standard Gaussian elimination 
in which a matrix is brought to reduced row echelon form rather merely 
to triangular form. In contrast to standard Gaussian elimination, entries 
above and below the diagonal have to be annihilated in the process of Gauss- 
Jordan elimination. It has been shown that the Gauss- Jordan elimination is 
considerably less efficient than Gaussian elimination with backsubstitution 
when solving a system of linear equations. Despite its higher cost, Gauss- 
Jordan elimination can be preferred in some situations. For instance, it 
may be implemented on parallel computers when solving systems of linear 
equations [2] . In addition, it is well suited for computing the matrix inverse. 

Applying Gauss- Jordan elimination to a given matrix A, we denote by A k 
the new matrix obtained after fcth step of Gauss- Jordan elimination. In the 
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present paper, we will show that each entry a^Ai 7^ j,j > k) in the matrix 
A k can always be expressed as a ratio of two determinants whose entries 
are from the original matrix A. In 2002, Gong et al. [1] first established a 
generalized Cramer's rule, which can be applied to a problem in decentralized 
control systems. However, their method is restricted to deal with a class of 
particular systems of linear equations. In [5], Hugo Leiva has presented 
another generalization of Cramer's rule, but the given formula is somewhat 
complicated. Different from the two methods mentioned above, our approach 
can also be used to directly construct one solution of AX = b. From this 
point of view, our method can give a generalized Cramer's rule whose form 
is completely different from the existing results. We also hope that it is useful 
not only as a theoretical tool, but also as a practical calculation methods in 
the linear algebra community. 



2. Main results 

Lemma 2.1. [3] If M is a square matrix and a, b, c, d are scalars, then 
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Before presenting the main result, we first offer a recursive description of 
Bareiss's standard fraction free Gaussian elimination [4]. 
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In what follows, in order to simplify the discussion, we also assume that 
the leading principal minors of a n x m matrix A are nonzero. 
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Theorem 2.2. Let A = (a^) be a n x m matrix with entries from an arbi- 
trary commutative ring and A k (0 < k < n) is defined as above. Bring A to 
reduced row echelon form by Gauss-Jordan elimination. Then after the kth 
elimination step, each entry a h Ai ^ j,j > k) in A k can be expressed as a 
ratio of two determinants whose entries are from the original matrix A. 

Proof. Consider the following three cases: 
1). Case 1: i > k, j > k. We shall show that 
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By (1), it is easy to see that the conclusion is true. To see this, let us use 
induction on the elimination step k as follows. 

(i) When k = 1, it is clear that the equality (2) holds. 

(ii) Now assume that the equality (2) is true for k. Then, when the elim- 
ination step is k + 1, we have 
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This proves the equality (2). 

2). Case 2: i = k,j > k. We shall claim that the below formula is true. 
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It is easy to prove this, since we want a k k 
elimination. 

3). Case 3: i < k,j > k. First, Let us construct the following determinant: 
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Next, we will claim that the following two recursion formulae hold. 
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Case 3-1. When % < k — 2, we have 
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Case 3-2. When i — k — 1, it follows that 
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The proof of the equality (6) : Since the row index of each element in the 
right-hand side of (6) is bigger than its column index, the formula (1) is still 
available. By (1), we get 
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Partition the above determinants into 4 submatrices respectively, as follows: 
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The last equality can be guaranteed by (4). 

A similar but somewhat more complicated method can be used to estab- 
lish the proof of (5). According to (1) and (4), we have 
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Afterwards, expand fl^-iife-i along the zth column, it follows that 
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Let M s be a square matrix whose determinant is M s . Since the minor M s 
obtained by expanding o^-i k-x a l° n g the zth column is exactly a minor of 
<4j , then one always can apply elementary row operations to a[j , such that 
the top left corner of a\j is exactly M s . Here, a[ k j = —a. 



According to Lemma 2.1, it follows that 
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And then, expand the above determinants along the ith column, we have 
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Now, we consider the third case: when i < k, j > k, the below equality 
holds. 
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(ii-3) When % — k, it follows that 
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Thus, the equality (10) holds. This completes the proof of Theorem 2.2. 
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According to the above results, we know that after kt\i Gauss- Jordan 
elimination step, each a^-(i 7^ j, j > k) in A k can be represented as a ratio of 
two determinants, as follows: 
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also easily construct one solution of AX = b. Thus, this method gives a 
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3. Conclusions 

As far as we know, the presented approach has not been published. This 
new method due to its distinct features can be used in a wide range of sci- 
entific and engineering problems. For example, it provides a feasible method 
to solve a system of linear equations with parametric coefficients by poly- 
nomial interpolation technique [8]. In addition, this method can be further 
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developed to give an explicit expression for the elements of the solution of a 
constrained linear systems of equations [7]. Finally, it can also be applied to 
solve some integer programming problems. 
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